Orbital order, anisotropic spin couplings, and the spin-wave spectrum of the 

ferromagnetic Mott insulator YTi0 3 
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Using a point-charge calculation of the electrostatic crystal field, we determine the non-degenerate 
orbital ground state of the ferromagnetic Mott insulator YTiOa, which is found to agree perfectly 
with experiment. Based on the orbital order, we obtain by perturbation theory an effective spin 
Hamiltonian that describes the magnetic superexchange between nearest-neighbor Ti ions. The su- 
perexchange Hamiltonian includes, in addition to the isotropic Heisenberg coupling, antisymmetric 
(Dzyaloshinskii-Moriya) and symmetric anisotropy terms, caused by the spin-orbit interaction on 
the Ti ions. We find ferromagnetic Heisenberg couplings for Ti-Ti bonds in the crystallographic 
ah planes, but antiferromagnetic ones for Ti-Ti bonds between planes, in contradiction with ex- 
periment (which gives ferromagnetic couplings for both). Difficulties in calculating realistic values 
for the isotropic couplings of YTiC>3 have been already reported in the literature. We discuss 
possible origins for these discrepancies. However, the much smaller values we obtain for the sym- 
metric and antisymmetric anisotropies may be expected to be reliable. We therefore combine the 
experimentally-deduced isotropic coupling with the calculated anisotropic ones to determine the 
magnetic order of the Ti ions, which is found to be in satisfactory agreement with experiment. 
Based on this magnetic order, we derive the spin-wave spectrum. We find an acoustic branch with 
a very small zone-center gap and three optical spin-wave modes with sizeable zone-center gaps. The 
acoustic branch reproduces the one reported in experiment, and the optical ones are in a satisfactory 
agreement with experiment, upon a proper folding of the magnetic Brillouin zone. 

PACS numbers: 71.10.-w, 71.27. +a, 75.10.Dg, 75.25.-l-z, 75.30.Dg 



I. INTRODUCTION 

The perovskite Ti oxides have attracted much interest 
since these strongly-correlated electronic systems possess 
orbital and magnetic degrees of freedom which are cou- 
pled together (for a review, see Ref. pj). A prominent 
member of this family is the ferromagnetic Mott insula- 
tor YTiC>3. The Curie temperature of this compound is 
Tq = 30 K and the ordered ferromagnetic moment, which 
is oriented along the crystallographic c axis, is 0.84 /l*b @- 
Another experiment |3J reported Tq — 27 K. Due to spin 
canting, there are also a small G-type antiferromagnetic 
moment along the a axis and a small A-type one along 
the b axis, which at T — 10 K amount to 0.08 {1b and 
0.05 hb, respectively (the c-axis ferromagnetic moment 
being 0.54 /xb at that temperature, which is extrapolated 
to 0.72 jiB at zero temperature) |3|. 

Several previous calculations aiming to explain YTiC>3 
and the doped series Lai-^Y^TiOs, respectively, have 
failed to achieve a consistent description of the ex- 
perimentally observed orbital and magnetic ordering 
EH A recent GGA+U ( generalized gradient ap- 
proximation + local Coulomb repulsion) study |7j has 
produced the correct orbital and magnetic ground state 
of YTi03, but has not provided quantitative estimates for 
the superexchange couplings between nearest-neighbor 
Ti ions. These couplings are required in order to un- 
derstand quantitatively the magnetic structure and the 
spin- wave spectrum observed in experiment 

From the microscopic point of view, YTi03 is quite 



similar to the antiferromagnetic Mott insulator LaTiOs: 
In both compounds there is a single electron in the 3d 
shell of Ti, and both have the same space group, Pbnm. 
It is worth noting in this connection that the magnetic 
structure of LaTi03 is also quite complicated: Experi- 
ment has indicated a predominant antiferromagnetic G- 
type moment along the a axis and a small ferromagnetic 
moment along the c axis and a recent theory 
has predicted a small A-type moment along the b axis. 
We have recently presented a detailed model for LaTi03 
, which proved successful in describing the orbital and 
magnetic ordering of that material and provided the su- 
perexchange couplings and the spin-wave dispersion mea- 
sured in experiment |lfjj |. Because of the apparent sim- 
ilarity between YTi03 and LaTi03, one may hope that 
the same model will explain the former as well. In this 
paper we carry out such investigation. Unfortunately, our 
model does not yield the correct isotropic Heisenberg su- 
perexchange coupling between nearest-neighbor Ti ions. 
The reason is that there are both ferromagnetic and anti- 
ferromagnetic contributions to that coupling, which have 
roughly the same order of magnitude. Our approxima- 
tions cannot resolve the competition between them to a 
sufficient precision. A similar problem has been reported 
in Ref. 0. 

Superexchange couplings are customarily derived per- 
turbatively, assuming that the hopping matrix elements 
are smaller than the on-site excitation energies. It 
has been found in other systems that such calcu- 
lations yield inaccurate values for the leading Heisen- 
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berg (isotropic) couplings, but are quite reliable for the 
much smaller anisotropic ones. It seems therefore rea- 
sonable to combine the experimental information on the 
isotropic couplings together with the calculated values of 
the anisotropic ones, in order to determine the magnetic 
structure of the ground state and the spin excitations. 
This procedure will be adopted in this paper. 

We begin in Sec. [H] with a point-charge calculation 
of the electrostatic crystal field due to all ions of the 
solid. The ground state of this crystal field determines 
the orbital order of the Ti ions. This orbital ordering 
agrees extremely well with the one detected experimen- 
tally 0, ^| . We then use a Slater-Koster parametriza- 
tion to compute the effective Ti-Ti hopping matrix el- 
ements. The Coulomb correlations on the doubly oc- 
cupied d shells are fully taken into account in terms of 
Slater integrals. Having thus obtained an effective mi- 
croscopic Hamiltonian for the Ti ions, we derive in Sec. 
IIHI the superexchange spin couplings between nearest- 
neighbor Ti ions, employing perturbation theory to lead- 
ing (second) order in the hopping matrix elements, and 
up to second order in the spin-orbit interaction on the 
Ti's. In this way we obtain, beside the isotropic su- 
perexchange coupling alluded to above, the antisym- 
metric (Dzyaloshinskii-Moriya) and symmetric superex- 
change anisotropics. Replacing the isotropic coupling 
by the experimentally-determined one, we calculate in 
Sec. IIVI the magnetic order of the classical ground state 
of the Ti ions and in Sec. [V] the spin- wave spectrum. 
The magnetic structure of the classical ground state is 
shown to be in satisfactory agreement with experiment 
The spin- wave calculation reproduces an acoustic 
branch of the spin-wave dispersion which has been de- 
tected by neutron scattering This branch has a very 
small zone-center gap, and is almost isotropic in the mag- 
netic Brillouin zone. In addition, we identify three op- 
tical spin-wave branches with considerable zone-center 
gaps. The experimental dispersion has been plotted as 
a single branch over the magnetic Brillouin zone (MBZ) 
of a pure ferromagnet However, YTiC>3 is a canted 
ferromagnet, for which the MBZ is four times smaller. 
We therefore prefer to re-plot the experimental spin- wave 
data according to the actual MBZ, i. e., to fold back the 
experimental data from the MBZ of the purely ferromag- 
netic case. When this procedure is adopted, one obtains 
a satisfactory agreement between the optical branches 
and experiment. In Sec. IVII we summarize our results 
and compare our picture of YTiOs with the ones given 
previously in the literature 0, 0, U E3 • 



II. THE MODEL 

A. The crystal field 

As is mentioned above, there is a single electron in the 
3d shell of the Ti ions in the ground state (YTi03 has the 
valences Y 3+ Ti 3+ (0 2 ~) 3 ). The unit cell, shown in Fig. [I] 
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FIG. 1: The crystallographic structure of YTiOa. The ten 
Ti ions, which constitute the twelve inequivalent nearest- 
neighbor Ti-Ti bonds are enumerated. For simplicity, oxygen 
octahedra are only shown around four Ti sites. Y ions from 
two layers are shown as small spheres. 



TABLE I: The structural parameters of YTi0 3 at T = 2 K 

E3- 



a 


5.32260 A 




0.f2133 


b 


5.695f 7 A 


yoi 


0.45702 


c 


7.59622 A 


X02 


0.690f0 


Xy 


0.97762 


yo2 


0.309f9 


VY 


0.07398 


ZQ2 


0.05770 



contains four Ti ions and twelve inequivalent nearest- 
neighbor Ti-Ti bonds. The crystal has the symmetry 
of the space group Pbnm (No. 62 in Ref. [Tj]). The 
structural data (taken at T = 2 K) are given in Tabled 
In order to use the symmetries of the space group, 
it is convenient to employ the orthorhombic orthonormal 
basis for the Ti-<i orbitals, 

\xy), |2z 2 ), \yz), \xz), \x 2 - y 2 ), (1) 

where the x, y and z axes correspond to the crystallo- 
graphic a, b and c axes. 

Using the structural data listed in Table HI we have 
calculated the spectrum and the eigenstates of Ti ion 
No. I (see Fig. employing a point-charge calculation 
of the static crystal-field Hamiltonian. This calculation 
uses the full Madelung sum over the crystal, which is 
evaluated as an Ewald sum |l6j |. It requires the second 
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TABLE II: The static crystal field for Ti 3+ (site 1): Spectrum 
and eigenstates in the orthorhombic basis for the d basis, see 



Eq. Q. 


-0.458 eV 


(-0.181, 


0.295, 


0.488, 


-0.542, 


0.590) 


-0.308 eV 


(-0.081, 


-0.412, 


0.529, 


0.653, 


0.343) 


-0.181 eV 


( 0.444, 


0.266, 


0.654, 


-0.017, 


-0.552) 


0.407 eV 


( 0.761, 


0.302, 


-0.231, 


0.222, 


0.477) 


0.540 eV 


(-0.430, 


0.762, 


-0.039, 


0.480, 


-0.040) 




FIG. 2: The orbital order of the Ti ions in YTi0 3 . 



moment, (V 2 ), and the fourth moment, (V 4 ), of the effec- 
tive ionic radius of the Ti 3+ -ion. We have used the values 
(r 2 ) = 0.530 A 2 and (r 4 ) = 0.554 A 4 [l7|. The results 
of the crystal-field calculation, which are listed in Table 
ITT1 exhibit a t% g splitting scheme where a non-degenerate 
ground state is clearly separated from the excited states. 
This ground state orbital, which gives rise to orbital or- 
dering, is given by the first line in Table 1TT1 and depicted 
in Fig. [21 It agrees very well with data obtained from nu- 
clear magnetic resonance (NMR) and polarized neutron 
diffraction experiments 



B. The Hamiltonian 

We next construct the microscopic Hamiltonian per- 
taining to the Ti ions, from which we obtain perturba- 
tively the superexchange spin couplings. The calcula- 



tion is carried out for a two-site cluster, consisting of two 
nearest-neighbor Ti ions, denoted by m and n. 

The unperturbed part of the Hamiltonian contains the 
static crystal field, H^ n , and the intra-ionic Coulomb 
correlations of a doubly occupied d shell, FL^ nn , 

TTO TTCf I TTC /o\ 

ran ran ' ran' \ I 

The perturbation calculation carried out below involves 
two sectors, which together span the Hilbert space of 
the cluster: A Ti 3+ sector, in which both Ti ions are 
trivalent, and a Ti 2+ sector, where one of the Ti ions is 
divalent (two d electrons on the same site) and the other 
is tetravalent (an empty d shell). The ground state of 
H^ nn belongs to the Ti 3+ sector, where both Ti ions are 
in the one-particle ground state of (modulo spin 

up or down on each site), leading to a four- fold degen- 
eracy of the ground state of the cluster. The spectrum 
of H^ nn is found by applying H^ n on the Ti 3+ sector, 
and both H^ ln and H^ n on the Ti 2+ sector. H^ n is 
parametrized in terms of the Slater integrals F 2 and F 4 
and the effective Ti-Ti charge-transfer energy U e e- 
This energy is the difference between the four-fold de- 
generate ground state of the cluster (which is the lowest 
level of the Ti 3+ sector) and the lowest level of the Ti 2+ 
sector (where H^ n and H^ nn are diagonalized simulta- 
neously). We use F 2 = 8F 4 /5 = 8.3 eV from an atomic 
Hartree-Fock calculation [13, and U e s = 3.5 eV from the 
analysis of the photoemission spectra and first-principles 
calculations [2fJ. We note that the charge-transfer en- 
ergy U c ff might have a considerable uncertainty, and in 
particular may be lower than 3.5 eV. For example, there 
is a strong resonance in Raman spectroscopy )2ll | at a 
laser frequency of 2.54 eV. Since experiment indicates 
that the resonance is mainly caused by processes involv- 
ing two Ti sites, it may well be that it yields a lower 
value for U c g. Further first-principles calculations and a 
comparison with optical-conductivity data are required 
in order to determine more precisely U e ff. 

The perturbation part of the Hamiltonian, V mn , con- 
sists of an effective Ti-Ti tunnelling term, i?j£!n, and the 
on-site spin-orbit interaction, H^ n , 

'mn mn ' van' V*-*/ 

The tunnelling Hamiltonian is given in terms of an effec- 
tive hopping matrix, t mn , between the m and the n Ti 
ions, 

H mn = ^ ] ^ ] ^mn^mia^nja ^" H. C., (4) 
ij a 

where d} mia ( d mia ) creates (destroys) an electron with 
spin a in the z-th eigen-orbital of H^ n at site number m 
(see Table ITT|l . The spin-orbit coupling is given by 

H% n = \J2h-Sk, (5) 
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where denotes the angular momentum operator of the 
Ti ion at the k site, S/. is its spin operator, and A is the 
spin-orbit coupling strength. We use A = 18meV [22$ . 

The dominant hopping process between two nearest- 
neighbor Ti ions is mediated via the oxygen ion which is 
nearest to both of them. Let tj£ be the hopping matrix 
element of an electron in the p orbital a on the oxygen ion 
into the i state of the Ti ion located at m. The effective 
hopping between the Ti ions is then given by 



fij — _ \ ^ fiat f jet _ ,ji 
ran \ / . m n ni 



(0) 



Here, A e ff is the charge-transfer energy, which is required 
to put an electron from an O ion on a Ti ion, and a 
denotes one of the three p orbitals on the oxygen (in 
orthorhombic coordinates), 



\x), \y), \z). 



(7) 



(Modifications of this basis due to the crystal field are 
ignored, since the crystal field splitting is expected to be 
small compared to the Ti-0 charge-transfer energy.) 

Using the structural data from Table [I] together with 
elementary geometric considerations, the Ti-0 hopping 
matrix elements can be expressed in terms of the Slater- 
Koster parameters V v da and V p d-n [23j . We use the values 
VU, = -2.3 eV, V pd7I = 1.1 eV, and A eff = 5.5 eV 
[20| . in conjunction with Eq. 10 to compute the effec- 
tive hopping matrices pertaining to the unit cell. The 
results are listed in Table ITTT1 which also gives the sym- 
metry properties of the hopping matrices between differ- 
ent Ti-Ti bonds. The four Ti sites of the unit cell form 
twelve nearest-neighbor Ti-Ti bonds which are inequiva- 
lent, i. e., they do not evolve from each other by Bravais 
translations. These bonds connect the ten Ti ions in- 
dicated in Fig. 2] By the symmetry operations of the 
space group Pbnm, the eight effective hopping matrices 
between Ti ions belonging to the same ab plane and the 
four matrices for Ti-Ti bonds along the c direction, re- 
spectively, can be expressed by a single matrix each. For 
example, all twelve hopping matrices are given by the 
two matrices for the Ti-Ti bonds ran = 12 (planar) and 
mn = 13 (inter-planar), respectively. 



C. The Ti-O hybridization 

Our model does not include the covalent contribution 
to the crystal field, arising from hybridization between 
the Ti-3d and 0-2p states. This mechanism mixes ex- 
cited states of the static crystal-field into the Ti 3+ ground 
state, i.e., there is an admixture of Ti 2+ states accom- 
panied by an admixture of holes on the oxygen sites. 

Following Refs. [2(J and |24|, we may estimate the ef- 
fect of the pd hybridization. When that hybridization 
is absent, the effective parameter U c s defines the en- 
ergy difference between the ground state of the Ti 3+ 



TABLE III: The effective Ti-Ti hopping matrices for the d 
eigen-orbitals of the crystal field from Table |H] values are 
given in eV. The rows and the columns are ordered beginning 
with the ground state of the crystal field (index 0), continuing 
with the first excited state (index 1), etc. The matrix tn is 
symmetric because of a mirror plane, see Ref. 

Planar 



tl2 


— tie — *25 


— *65 


— *34 — 


^38 — *47 


— *87 




" -0.062 - 


0.206 


0.033 


-0.026 


-0.012 




0.007 - 


0.015 


0.006 


-0.086 


0.114 




0.130 - 


0.077 


-0.125 


0.149 


-0.203 




-0.202 


0.030 


0.092 


0.453 


-0.632 




-0.036 


0.008 


0.024 


0.031 


-0.044 



Inter-planar 



— ^24 — *39 


— '410 








0.086 - 


0.009 


0.101 


-0.024 


-0.085 


-0.009 


0.160 


0.043 


0.126 


0.227 


0.101 


0.043 


0.119 


-0.048 


-0.159 


-0.024 


0.126 


-0.048 


-0.107 


-0.263 


-0.085 


0.227 


-0.159 


-0.263 


-0.607 



sector and the lowest state of the Ti 2+ sector in a two- 
site cluster consisting of two Ti ions. When the pd hy- 
bridization is present, these two types of d states corre- 
spond to two bands, from which two pd hybridized bands 
evolve according to the covalent crystal field. These hy- 
bridized bands have, in general, significant dispersion: 
Their peak-to-peak separation is given by the band gap 



-Egap=L8 eV |2fJj, and the distance between the band 
edges is given by the optical gap £ , O pt=l-0 eV which 
is experimentally observed as the Mott gap. The mean 
bandwidth between the two pd hybridized bands is then 
W = E gap —E opt = 0.8 eV. These bands are not as disper- 
sive as in the case of LaTiOs, where the mean bandwidth 
is W = 1.4 eV H|. 

Nevertheless, given this dispersion of the bands one 
may wonder whether a localized picture is appropriate, 
even approximately, for the YTiC>3 system. In order to 
study this point, we have analyzed the covalent crystal 
field of a cluster consisting of a single Ti ion, and the six 
oxygen ions predominantly hybridized with it (the calcu- 
lation has been carried out for Ti number 1 in Fig. ^| . 
This is accomplished by diagonalizing the Hamiltonian 

cf 



h c + m 



(8) 



for a TiOe-cluster. Here H describes the static crystal 
field, H c is the Coulomb interaction, and 
pd-tunnclling, 



is the 



rrtun 

11 pd — 



nicta 



t ia d^ n 



H.c, 



(9) 



where p na a destroys an electron on the n-th oxygen site 
with spin a in the a-orbital, given in Eq. 0. As in 
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TABLE IV: The combined static and covalent crystal field for 
Ti 3+ (site 1): Spectrum and eigenstates in the orthorhombic 
basis for the d basis, see Eq. Q. 



-0.673 eV 


( 0.187, 


-0.340, 


-0.438, 


0.583, 


-0.564) 


-0.519 eV 


(-0.028, 


-0.350, 


0.573, 


0.622, 


0.402) 


-0.409 eV 


( 0.459, 


0.274, 


0.634, 


-0.050, 


-0.557) 


0.737 eV 


( 0.751, 


0.342, 


-0.280, 


0.188, 


0.453) 


0.865 eV 


(-0.435, 


0.755, 


-0.036, 


0.485, 


-0.072) 



the calculation of the Ti-Ti hopping amplitudes, the pd 
hopping amplitudes, t|", are expressed in terms of the 
Slater-Koster parameters V v da and V p dn, using the struc- 
tural data of Ref. . 

The entire space of the basis states of the TiOg-cluster 
consists of a Ti 3+ sector where the p orbitals are all oc- 
cupied, and a Ti 2+ sector where there is a hole in one 
of the p orbitals. The eigenstates of the Hamiltonian (JSJ 
have the form 

=\/2 — ^jd^ + V^TIrf 2 ), (10) 

where rid is the occupation number of the Ti-d shell (1 < 
nd < 2), jd 1 ) is a state with a single electron in the 
d shell and fully occupied p shells on the surrounding 
oxygen ions, and |<i 2 ) is a state with two electrons in the 
d shell and a hole in the p shell of one of the oxygen 
ions. We find that in the ground state rid — 1.330 , i. e., 
there is a p hole on one of the neighboring oxygens with 
probability of 33.0 %. 

This calculation allows for the analysis of the eigen- 
states of the combined static and covalent crystal fields. 
Projecting the five lowest eigenstates onto the Ti 3+ sec- 
tor (which corresponds to the states jd 1 )), gives to a very 
good approximation the same eigenstates as for the static 
crystal field alone, as can be seen by comparing Table llVl 
with Table [U] This finding explains why, despite the 
admixture of Ti 2+ states \d 2 ), the agreement with the 
experiments of Refs. 01 an d remains perfect. In- 
deed, these experiments measure the Ti 3+ part, jd 1 ), of 
the combined static and covalent crystal field, and ap- 
parently are not sensitive to the Ti 2+ admixture \d?)- 
Table Hvl also shows that the t2 g splitting remains almost 
the same as in the absence of the covalent contribution, 
whereas the distance between the t^g and e g energies is 
enhanced. 

We now discuss the crystal-field gap and the t2 g split- 
ting scheme as obtained from our calculation and from an 
alternative calculation 123. in relation with an analysis of 
the optical conductivity |2fij and Raman data [2ll |. Our 
static crystal-field calculation yields a non-degenerate or- 
bital ground state separated by «0.15 eV from the first 
excited state and a second excited state separated by 
«0.13 cV from the first excited one (see Table HTJ. This 
t2 g splitting scheme results from the orthorhombic distor- 
tion of the crystal and from the distortion of the oxygen 



octahedra. We have estimated that the covalent crystal 
field reduces the gap between the first and the second 
excited states (to about 0.11 eV according to Table llVjl. 
while the gap between the ground state and the first ex- 
cited state remains practically the same. A more precise 
calculation of the covalent crystal field |2^|, which takes 
into account two additional effects, the pp hybridization 
and the Ti 1+ admixture [27]], gives «0.19 eV for the gap 
between the ground state and the first excited state and 
«0.14 eV for the gap between the first and the second 
excited states (the Ti 1+ admixture means that there is 
also an admixture of d 3 states to the ground state). This 
result is in better agreement with the data of optical con- 
ductivity |2^| and Raman spectroscopy j2l|, which show 
that the first orbital excitation is centered around 0.2- 
0.25 eV. 

Since it is extremely complicated to include in the mag- 
netic superexchange calculation the hopping between the 
pd hybridized states, our calculations below contain only 
the hopping between the Ti 3+ states. The results listed 
in Tablc llVl which show that the projections of the eigen- 
states of the combined static and covalent crystal fields 
onto the Ti 3+ sector are almost the same as in the static- 
only case, ensure that the Ti 3+ ground states we use are 
an appropriate starting point for the superexchange cal- 
culation. 



D. The magnetic moment 

The calculation of the magnetic structure detailed be- 
low yields the directions of the magnetic moments in the 
ground state, but does not determine the magnitude of 
the moment. However, one can estimate that magni- 
tude by diagonalizing together H^ n and H^ n for a sin- 
gle Ti 3+ ion. The eigenstates of this combined Hamil- 
tonian are symmetric or antisymmetric with respect to 
time-reversal, leading to five Kramers doublets for the 
single Ti 3+ ion. We use those doublets to find the ex- 
pectation values of the angular momentum. By choosing 
the largest possible polarization of the magnetic moment 
along the z axis (that direction is the leading one of the 
observed moment Q) out of all the linear combinations of 
the ground-state doublet, we find (^+2s|) /i B = 0.91 /z B . 
This partially explains the reduction of the observed or- 
dered moment as compared to 1 hb- 

The Ti-0 hybridization hardly affects the magnetic 
moment. For the parameters used here, the admixture 
of spin and spin 1 Ti 2+ states into the ground state of 
the covalent crystal field reduces the ordered moment by 
only w 0.1 %. 



III. THE SUPEREXCHANGE COUPLINGS 

Our aim is to obtain from the full Hamiltonian, H mn = 
Hf nn + V mn , an effective spin Hamiltonian, h mn , which 
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acts within the Hilbert space of the four-fold degenerate 
ground state of the unperturbed Hamiltonian H® nn . 

Since the Hamiltonian is invariant under time-reversal, 
there are no single-ion terms, and consequently the effec- 
tive spin Hamiltonian, to second-order in the spin vari- 
ables, takes the form 

hrnn — S m • A mn • S n , (H) 

where A mn ( = A l nm ) is the 3x3 superexchange ma- 
trix. This matrix may be decomposed into a symmetric 
part and an antisymmetric one. The three components of 
the latter constitute the Moriya vector D mn (= — D„ m ). 
Extracting further the isotropic part of A mn , i.e., the 
Heisenberg coupling J m „, the effective spin Hamiltonian 
is cast into the form 

h m n — ^mnSm'Sn+D rml -(S m XS Il )+S TO -^ ml -S Jl . (12) 

Here, A s mn represents the symmetric anisotropy. Due to 
the space-group symmetries, all three types of magnetic 
couplings belonging to the eight planar Ti-Ti bonds may 
be obtained from those of a single bond, and so is the 
case for the four inter-planar bonds, see Ref. @. 

The various magnetic couplings appearing in Eq. 1|12|) 
are obtained by perturbation theory to leading order in 
V mn , namely, to second order in the hopping t mn and to 
first and second order in the spin-orbit coupling (scaled 
by A). The formal expressions of the perturbation ex- 
pansion are documented in Ref. 0. The Heisenberg 
isotropic exchange [the first term in Eq. I|12fl ] is inde- 
pendent of A. A systematic description of the magnetic 
anisotropics due to the spin-orbit interaction requires 
both the first and the second order processes in A |29j . 
The technical reason being that the expectation value of 
the cross product in the second term of Eq. (|12[1 is, in fact, 
also of order A, so that altogether the Dzyaloshinskii- 
Moriya interaction contributes to the exchange energies 
in at least second order in the spin-orbit coupling. We 
neglect terms in which there appear two Ti 2+ interme- 
diate states in the perturbation expansions. These are 
smaller than the ones we keep, by an additional factor 
of ~ A cf /U cfi = 0.043, where A cf = 0.150 eV is the gap 
between the ground state of the single-particle crystal 
field and the first excited state, see Table [H] The de- 
tailed calculation of the various terms is lengthy, albeit 
straight-forward. More details are given in Ref. 91. The 
values we obtain, using the parameters cited above, are 
listed in Table E| 

As is seen from Table [3 the value we find for the in- 
plane Heisenberg coupling, Jj9, roughly agrees with the 
experimental one, «— 3 meV |3j. However, the calculated 
inter-plane coupling, J13, is positive, in contradiction 
with experiment (which yields for that coupling rj -3 
meV, too). Generally speaking, there are antiferromag- 
nctic and ferromagnetic contributions to the isotropic 
Heisenberg couplings. The former arise when the inter- 
mediate Ti 2+ states of the perturbation expansion are 
singlets, and the latter when they are triplets. Separat- 
ing these competing contributions, we find Jf 2 =21.481 



TABLE V: The calculated single-bond spin couplings (in 
meV). The Moriya vectors are given including the cor- 
rections T)' mn , which are of order A 2 . The symmetric 
anisotropies are given as A„ n = (4" n , A% n , A™ n ) and 
A-mn = {Amni A^ n , A^„) for the diagonal and off-diagonal 
entries, respectively. 

Heisenberg couplings 
J12 = -3.870, J13 = 2.772 
Moriya vectors 
D12 = (1.776, -0.938, -0.325), D13 = (-0.671, 0.189, 0) 
Symmetric anisotropies 
A? 2 = (0.175,-0.011,-0.160), A? 3 = (-0.145,-0.024, 0.010), 
A°2 = (0.044, -0.131, -0.313), A°% = (0, 0, -0.153) 



meV, J 12 =-25.351 meV, Jf 3 =12.008 meV, and J 13 =- 
9.237 meV, namely, the antiferromagnetic and the ferro- 
magnetic contributions are roughly the same, ft is worth 
noting that in the case of LaTiO*3 9], the contribution of 
the singlets dominated the one of the triplets, and indeed 
our calculation of that compound has yielded reliable val- 
ues. Unfortunately, in the case of YTi03 the balance 
between these competing contributions is too sensitive 
to be resolved by our model approximations. This del- 
icate balance in the case of YTiC>3 is also reflected in 
the overall rather small value of the total isotropic cou- 
pling, ss-3 meV, (whereas it is 15.5 meV in the case of 
LaTiOs). This value is also very sensitive to the parame- 
ters used. For example, taking U e g « 1.6 eV would have 
changed the sign of J13. In contrast, a sign change in the 
case of LaTiC>3 requires the considerably lower value of 
U e ff w 0.6 eV. We discuss in Sec.|^various difficulties en- 
countered in obtaining realistic values for the Heisenberg 
couplings of YTiC>3 which have been previously reported 
in the literature. 

Had we used the values listed in Table |Vl we would 
have obtained a predominant A-type antiferromagnetic 
order for YTiC>3, which sharply contradicts the experi- 
ment Q . However, the fact that our values for the leading 
(isotropic) superexchange couplings do not agree with ex- 
periment does not necessarily mean that the anisotropic 
ones are not reliable. In the case of the cuprates, for 
example, it has been found (by comparing with ex- 
act diagonalizations) that while the isotropic couplings 
calculated by perturbation theory were inaccurate, the 
anisotropic ones were accurate enough |ll| . Since the 
latter determine the directions of the spins in the classi- 
cal ground state and the spin-wave gaps, a way to test 
our anisotropic superexchange couplings is to examine 
those properties. In order to do so, we replace in the 
following the isotropic couplings by the experimentally- 
deduced ones, J12 = J13 = -2.75 meV @, while using 
for the anisotropic couplings the values given in Table [V] 
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IV. THE MAGNETIC STRUCTURE 

The single-bond spin Hamiltonian, Eq. I|12f> . is the ba- 
sis for the magnetic Hamiltonian, from which the mag- 
netic order of the classical ground state follows. To con- 
struct the latter, the entire Ti-lattice is decomposed into 
four sublattices. The four sublattices are enumerated ac- 
cording to the numbers of the four Ti ions per unit cell 
shown in Fig. ^ (sublattice i = 1 corresponds to Ti ion 
1 and its Bravais translations, etc.). Assigning a fixed 
magnetization (per site) to all the spins within each sub- 
lattice, Mi , one sums over all bonds which couple the four 
sublattices, to obtain the macroscopic magnetic Hamil- 
tonian in the form 

//>• = E [ J « M < ■ x l + D § ■ ( M * * M i) + AI < ' r -.. ' M j] > 

(13) 

where ij runs over the sublattice pairs 12, 13, 24, and 34 
of Fig. n This summation procedure gives rise to the 
macroscopic magnetic couplings: 1^ is the macroscopic 
isotropic coupling, D^- are the Dzyaloshinskii vectors 
which are the macroscopic antisymmetric anisotropies, 
and r f • are the macroscopic symmetric anisotropy ten- 
sors. The relations between those macroscopic cou- 
plings and the microscopic single-bond couplings and 
the inter-relations among the macroscopic couplings, dic- 
tated by the symmetries of the space-group, are discussed 
in Ref. 0. 

We now minimize Hm, and find the various sublat- 
tice magnetizations. We assume that all four vectors Mj 
have the same magnitude, M, but differ in their direc- 
tions. Since Eq. (|13fl is quadratic in M, the minimization 
will only yield the directions of these vectors, and not the 
value of M . To simplify the procedure we use group the- 
ory. According to the space group Pbnm symmetries, 
there are four possibilities for the symmetry of sublattice 
magnetizations of the classical ground state, as listed in 
Ref. |3(|. Having checked each of them, we have con- 
cluded that only one of these possibilities has the lowest 
energy. We then find that the classical magnetic ground 
state has the following structure: The x components of 
the magnetizations order antiferromagnetically, in a G- 
type structure (where the four sublattices actually reduce 
to two). The y components order antiferromagnetically 
as well, but in an A-type structure. Finally, the z com- 
ponents of the magnetizations order ferromagnetically. 



This structure agrees with experiment |3j . This is a non- 
trivial result caused by the anisotropic spin couplings. 
Given only the ferromagnetic Heisenberg couplings, the 
ferromagnetic moment could also be oriented along the 
x or the y axis, see Ref. |3fJ. 

The magnetic structure can be specified by expressing 
the four magnetizations in terms of two canting angles, 
tp and z?, see Tables IvTI and IVTil The angle measured 
with respect to the z axis, is proportional to the spin- 
orbit parameter A (as found by varying this parameter), 
while the angle tp is almost independent of it. Indeed, 

TABLE VI: The structure of the magnetic order in YTi0 3 , 
characterized in terms of the sublattice magnetizations Mi in 
the classical ground state (normalized to the value M), which 
are expressed by the canting angles tp and 



x components: G-type 


-Mf = 


Mf = Ml = -Mf = Mcos<£sintf 


y components: A-type 


-Mf = 


-Ml = Ml = Ml = M sin tp sin ■& 


z 


components: ferromagnetic 


M{ 


= M| = Mf = Ml = Mcostf 



one may verify that for an infinitesimally small A, there 
is just a ferromagnetic order along the z axis. As A in- 
creases, so does d, causing an increasing canting of the 
magnetizations. However, the projection of the magnetic 
moment onto the xy planes remains almost perpendicu- 
lar to the rotation axis of the magnetization, and hence tp 
is practically unaffected by the value of A. Interestingly 
enough, the magnetic structure of LaTi03 can also be 
described in terms of such canting angles. However, in 
that case both (p and $ (the latter measured with respect 
to the xy planes) are proportional to A [9|, leading to a 
(mainly) G-type order along the x direction which would 
have occurred even for an infinitesimally small A. 

It is worth noting that using naively the procedure out- 
lined above to obtain the energy of the classical magnetic 
ground-state might yield non-systematic contributions 
up to fourth order in the spin-orbit coupling A |9tl29|. To 
exemplify this point, we consider the expectation value 
of Hm, expressed in terms of the canting angles tp and 
0, 
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(H M ) = [\°:] 2(/ 12 +/ 13 )cos 2 tf 

[A 2 : ] - 2(J 12 + Z 13 ) cos 2 (f sin 2 -d + 2(I 12 - 7 13 ) sin 2 tp sin 2 $ 



[A 4 : 



D® 3 V ) cos 95 cos tfsini? - sin cos t? sin ± 2(rf| + rff) cos 2 1? 



4I?P 2 Z cos tp sin sin 2 ^ — 2T\ 2 sin 95 cos 1} sin $ 
2 (rff + rff) cos 2 tp sin 2 ± 2(rf -r»^« 2 



r^|) sin tp sin — 4r i g cos 99 sin 99 sin z9. 



(14) 



TABLE VII: The macroscopic magnetic couplings in meV, 
the resulting canting angles of the magnetizations in the clas- 
sical ground state, and the resulting values of the ordered mo- 
ments (normalized to the value M). Three coefficients of the 
macroscopic symmetric anisotropics are taken into account 
(see text). 

Isotropic couplings 

J12 = -2.750, J13 = -1.375 

Dzyaloshinskii vectors 
D?2 = (0,-0.938, -0.367), D? 3 = (-0.335,0.094,0) 
Macroscopic symmetric anisotropics 

rff = -0.160, rff = 0.005, rf| = 0.044 

Canting angles 
ip= -44.17°, = 7.55° 
Ordered moments 
M = (±0.094, ±0.092, 0.991) M 



We take the contributions up to the order A 3 into ac- 
count, i. e., we exclude from the calculation of the classi- 
cal ground state the coefficients rff, rff, T\\, T^, Dj^, 
and the A 2 correction of Df 2 z ■ This procedure yields the 
macroscopic magnetic couplings listed in Table IVTT1 Us- 
ing these couplings we have calculated the canting angles 
tp and and the ordered magnetic moments. These re- 
sults are also listed in Table IVTT1 

In order to compare our magnetic structure with the 
one found experimentally, we normalize the moments to 
1 fiB- Then, according to Ref. |3J, experiment gives M = 
(±0.149, ±0.085, 0.985) /x B , with relative errors of 15 % 
for the G-type moment, 25 % for the A-type moment, 
and 2 % for the ferromagnetic moment. The calculated 
values are within a single error bar except for the G-type 
moment for which we obtain a value which is 37 % lower 
than the experimental one (i.e., within the 3<r range of the 
measurement). Thus, the calculated magnetic structure 
is in reasonable agreement with experiment. 



V. THE SPIN- WAVE SPECTRUM 

A. The spin-wave Hamiltonian 

The calculation of the spin-wave dispersion is carried 
out analogously to the case of LaTi03 [|| , since all sym- 
metries are the same for both YTi03 and LaTi03. As in 
the calculation of the classical magnetic ground state, we 
combine the experimental Heisenberg couplings J12 = 
J13 = -2.75 meV with our calculated anisotropic 
couplings which contribute systematically to the classi- 
cal ground-state energy [see Eq. (|14fl and the following 
discussion] . 

Since the classical magnetic ground state is character- 
ized by four sublattices, we will obtain four branches in 
the spin-wave dispersion. The first step in the standard 
calculation of spin- wave dispersions is the rotation of the 
local coordinates at each sublattice, i, such that the new z 
axis will point in the direction of the corresponding sub- 
lattice ground-state magnetization, Mi. This rotation 
still leaves the freedom to choose the new local x and y 
axes, i.e., to rotate the new coordinate system around 
its z axis. Denoting the new local coordinate system by 
x'i, y'i and z[ (i = 1,2,3,4), we find that the convenient 
choice for our purpose is 

y[ = — , x =y xi, M=Mi, 

M uii 11 

mi =yf(MV Y + (M? f. (15) 

In the rotated coordinate system the spin Hamiltonian, 
comprising all three types of magnetic couplings, takes 
the form 

h = S' m ■ A' mn ■ S' n , (16) 

(ran) 

where the primes denote the rotated quantities. In par- 
ticular, A' mn is the 3x3 superexchange matrix in rotated 
coordinates. 

Since we consider only the Ti ions, it is convenient to 
use a coordinate system in which the Ti ions occupy the 
sites of a simple cubic lattice, of unit lattice constant 
(this picture is the appropriate one for comparing with 
the experimental spin- wave data Q, as discussed in the 
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next subsection). It is also convenient to use a coordi- 
nate system in which nearest-neighbor Ti ions are located 
along the axes (namely, to rotate the orthorhombic coor- 
dinates by —45° around the z axis, see Fig. QJ. Accord- 



ingly, our magnetic unit cell is spanned by the vectors 
(1,1,0), (1,-1,0), and (0,0,2), and the corresponding 
magnetic Brillouin zone (MBZ) is defined by 



\qx + q y \<K, |<fe|<2" 

The resulting spin- wave dispersion (for more details of the derivation, see Ref. |9() consists of four branches, 



fii(q) = (Ci + C% cosq 2 ) 2 - \Cz\ 2 cos 2, q z + \c\ l^cos^ + cosg y ) 2 - |C 3 cosq x + C 3 '* cosq y 

-(cosq x + cos q v )W(- cosg 2 ), 
n 2 (q)=fi 2 (q + Q), with Q = (0,0,7r), 
ft|(q)=n 2 (q + Q'), with Q' = (tt.tt.O), 
ft|(q)=fi?(q + Q"), with Q" = Q + Q' = (7r,7r,7r), 



-.11 



-.11* 



where 



I II * 

^ 2 (co Sfe )=4[(C 1 - ^cosg 2 ) 2 - I^IWj,] [|C|| 2 - ( C 3+^3 ) 2 
{Ci*c\ + Cicf) cosq z + (Ci - Ci cos q z ){c\ +cf) 



(18) 



(19) 



Each of the spin- wave branches has tetragonal symmetry, i.e., Qi(q x , QyiQz) = ^i(%> 1x,<lz) 

= Q,i(—q x ,q y ,q z ) — Qi(q x , —q y , q z ) = Sli(q x ,q yi —q z ). The coefficients in Eqs. <|18ll and 119|1 are linear combinations 
of the coefficients C mn (£), 

C 1 =2C 13 (l) + 4C 12 (l) = Cr, Ci = 2C 13 (2) = C 2 X *, cj = 2C 12 (2), C x = 2C 13 (3), c| = 2C 12 (3). (20) 
These are given by combinations of the superexchange matrix elements (A' mn ) al3 , 

n {X\ = — L( a' Y z 

mn V / o v mnJ 5 



Cnm(2)- J [(^mn) 2 ^ + (4»n) W + '((^mn) 91 ~ (^mn)^)] ! 
Cmn(3)= 4 [(^-mrx)^ — {■^mn) VV + H^ron)^ + (^mn) 1 ')]- 



(21) 



The explicit expressions are not reproduced here since 
their expressions are very long. 

Equations (|18fl contain our final result for the spin- 
wave spectrum of YTi0 3 . Evidently, the details of the 
spectrum can be obtained only numerically: One has to 
write the spin- wave coefficients, Eqs. <|20[1 . in terms of 
those appearing in Eqs. (|21|l . and express the latter in 
terms of the original coefficients of the spin Hamiltonian 
(|12fl . These results are then used in constructing the 
dispersion. 

When the spin-orbit parameter A is set to zero the 
coefficients appearing in Eqs. (|18|) simplify to 

C 1 =— 2J 12 — J 13 , C2 — J 13 , cj = J 12 , C% — C 3 =0, 

(22) 

where Ji 2 < is the in-plane Heisenberg coupling, and 
J13 < is the Heisenberg coupling between planes. In 



that case 

fi?(q)=[2J 12 + J 13 - J 12 (cosg x + cosg 1/ ) - J 13 cos^] 2 , 

^l(q)= [2 J12 + <As - J 12 ( cos Qx + cos Qy) + J 13 cos Qz] > 

fi 2 (q)= [2 J 12 + J 13 + J 12 (cos q x + cos g v ) - J 13 cos gj 2 , 

nf(q)=[2J 12 + J 13 + J 12 (cos 9a; +cosg y ) + J 13 cos^] 2 . 

(23) 

Only fii(q) vanishes at the zone center and is hence 
termed the acoustic mode. The other branches have 
gaps at the zone center: fi 2 (0) = 2| Ji 3 |, ^3(0) = 4| Ji 2 |, 
514(0) = 4| J12I + 2 1 Ji 3 |, and are hence termed the optical 
modes. Indeed, when only the ferromagnetic couplings 
J12 and Ji 3 are kept (i.e., for A = 0), the magnetic unit 
cell contains only one Ti ion, corresponding a simple cu- 
bic lattice. The Brillouin zone is then four times as large 
as the Brillouin zone of Eq. (|17fl . By "folding out" the 
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three optical modes into the larger Brillouin zone, one 
reproduces the usual gapless dispersion of the pure (fer- 
romagnetic) Heisenberg model. At finite values of the 
spin-orbit coupling all modes have gaps at the zone cen- 
ter, but the one of f2i is much smaller than those of the 
other three modes. 



anisotropy parameter A = 0.8 meV (which expresses the 
symmetric anisotropies which are allowed in a cubic situ- 
ation), and a zone-center gap A = 0.093jjr = 0.02 meV, 

the authors of Ref. have fitted the measured neutron- 
scattering data numerically onto the dispersion, 



B. Numerical results for the spin- wave dispersion 

1. The acoustic branch 

Using an isotropic ferromagnetic nearest-neighbor 
Heisenberg coupling J = — 2. 75 me V for all bonds, an 

I 



fi(q) ~W | J|[3 — (cos q x +cosq y + cos q z )] + A + A(l — cosq x W | J| [3— (cos q x + cosq y +cosq z )} + A +A(l — cosq y ). (24) 



The numerically-fitted zone-center gap is extremely 
small. On the other hand, Ref. reports an upper 
bound, 0.3 meV, for the gap A. We find for the acoustic 
branch 



onal anisotropy of the acoustic branch is found to be 



Ai = (^(0) = 0.326 meV. 



(25) 



This value for the zone-center gap roughly agrees with 
the upper bound according to Ref. |3j. A more severe 
discrepancy concerns the anisotropy parameter A. This 
parameter implies that the diagonal and off-diagonal en- 
tries of the symmetric anisotropy tensors are given by 

A i2 = |(A A, 0) = (0.4, 0.4, 0) meV, 



L 13- 

. od 



=(0,0, A) = (0,0,0.8)meV, 



Ag=|(0, 0, -A) = (0, 0, -0.4) meV, 
A? 3 d =(0,0,0). 



(26) 



This result is in contrast with our calculated values for 
the symmetric anisotropy tensors according to Table Ivl 
However, as is noted above [see Eq. <jT4*|) ] and also 
elsewhere both antisymmetric as well as symmetric 
anisotropies contribute to the same order in the spin- 
orbit parameter to the classical ground-state energy and 
hence to the spin-wave gap and its dispersion. In other 
words, it is not possible to express all anisotropies in 
terms of a single parameter and the zone-center gap, or 
alternatively, it is not possible to deduce the strength of 
the Dzyaloshinskii-Moriya interaction directly from the 
spin-wave dispersion. A proper procedure is to compare 
the full measured and calculated dispersions. Figures 
(a)-(c) show such a comparison. The agreement of the 
calculated acoustic branch with the experimental func- 
tion of Eq. I|24|) is satisfying (though it is quantitatively 
not as good as in case of LaTiOa). The calculated tetrag- 



v ' ' 2> = 97 74 % 

fii(i,o,o) y ^ 4/0 - 



(27) 



This value is smaller compared to the one found for 
LaTiOa |9j, mainly because in the present case the 
Heisenberg couplings are taken to be isotropic over the 
lattice. 



2. The optical branches 

There is experimental evidence for optical spin-wave 
branches, though this has not been pointed out explic- 
itly in Ref. 0. There, a plot of the dispersion along 
the (0,0,1) direction (in the pseudocubic coordinates also 
used here) is shown. The plot range includes wave vec- 
tors (in our notation) q = 7r/2(0,0, 1)£ with £ == 0...2. 
One may note however, that for £ > 1 these wave vectors 
are located outside the (first) MBZ, see Eq. I|17fl . The 
reason is that Q = (0, 0, n) [as well as Q' = (tt, 7r, 0)] is 
a site on the reciprocal lattice of sublattice No. 1 and 
thus, is equivalent to zero wave vector. Were YTi03 a 
ferromagnet without any spin canting, then all four mag- 
netic sublattices would have combined to form a single 
lattice, and Q, Q', and Q" = Q + Q' would not be sites 
of the reciprocal lattice. In this hypothetical case the 
MBZ would be given by \q x \, \q y \, \q z \ < n. In the ac- 
tual case, with the spin canting, the MBZ is four times 
smaller than that. We prefer to fold back the data from 
Ref. Q to this smaller MBZ, i. e., not to consider only 
the experimental fit function f2(q), but also 0(q + Q), 
fi(q + Q'), and n(q + Q"). We have plotted Figs. El (a)- 
(c) using the back-folded MBZ. The experimental curves 
which are folded back into the MBZ of the canted ferro- 
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(a) 

c ,7/2:1,1.1* 



17.5 




0.2 0.4 0.6 0.S 1 



FIG. 3: The spin-wave dispersion along selected directions in 
the magnetic Brillouin zone. Panels (a)-(c) show the four 
branches Q»(q) of our calculated dispersion (solid curves), 
the single branch fi(q) (dashed curves) which has been fit- 
ted onto neutron-scattering experiments, Eq. (1241 . according 
to Ref. 0, and three branches, which are obtained from f2(q) 
by folding it back from the MBZ of the uncanted ferromag- 
net into the smaller MBZ of the canted ferromagnet (dot- 
ted curves, see text), (a) The dispersion along the direction 
(1,1,1); (b) the dispersion along (1,0,0); (c) the dispersion 
along (0,0,1). 

magnet are in a satisfying agreement with our calculated 
optical branches. In fact, the signature of the optical 
spin-wave mode with the largest zone-center gap can be 
seen in the plots of Ref. There, it is related to the 
wave vector Q" (which is equivalent to zero wave vector) 
and has an energy of about 18 meV. 

We suggest to re-analyze also the spin-wave data on 
LaTiC>3 [l(|, i- e., to fold the experimentally-deduced 
dispersion fi(q) from the MBZ of the antiferromagnet 
without spin canting back to the MBZ of the canted an- 
tiferromagnet (which is half as large), considering also 
17 (q + Q). Then one obtains the result that the optical 
branches which have been calculated in Ref. [j| are con- 
sistent with the experimental fit function of Ref. ^} f° r 
the dispersion. 

In the following, we summarize the properties of the 
calculated optical branches of YTiC>3 . They have consid- 
erable zone-center gaps, 

A 2 =17 2 (0) = 5.815 meV, A 3 = 17 3 (0) = 11.721 meV, 
A 4 =17 4 (0) = 17.214 meV. (28) 



(b) 

q=«(1 ,<!,□)<; 



17.5 




□ .2 0.4 0.6 0.8 1 

i 



(c) 

q=jr/2(0 ; 0,1)f 



17.5 



Ofq) 

[meV] 




□ .2 OA 0.6 0.S 1 



Two of the calculated optical branches have considerable 
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tetragonal anisotropies, 



n 2 (o,o,f) 
o 2 (f,o,o) 

»4(0,0,f) 

n 4 (f,o,o)' 



=35.48%, 



=100.77%. 



" 3 (0, 0,f) 
03(5,0,0) 



= 163.29%, 



(29) 



As is seen in Figs. 02 (a)-(c), all four spin-wave modes 
are highly non-degenerate over a wide range of the MBZ. 
In contrast, in the case of LaTiC>3 we have found |j that 
the four modes constitute two pairs of quasi-degenerate 
branches. The reason for this difference between the 
two systems is related to the smallness of the angle ip 
in LaTi03 (as opposed to its significant value in the case 
of YTi0 3 , see Table lyTT)) . In LaTi0 3 there is a nearly full 
translational symmetry, leading to the quasi-degeneracy 
of the modes. 



VI. SUMMARY AND DISCUSSION 

We have presented a detailed model that aims to de- 
scribe the orbital and the magnetic orders in YTiC>3. 
While the orbital order that we have calculated turns out 
to agree very well with experiment, this is not the case for 
the magnetic superexchange couplings: To the lowest or- 
der in perturbation theory, we find that the approximate 
isotropic coupling between the ab planes is antiferromag- 
netic, while experiment indicates that it is ferromagnetic. 
This discrepancy is apparently due to a strong competi- 
tion between ferromagnetic and antiferromagnetic contri- 
butions to that coupling. The approximations we employ 
are not sensitive enough to resolve successfully this com- 
petition. In particular, the neglect of exchange processes 
which involve double p holes in the intermediate states, 
and of the Ti 2+ admixture into the ground state caused 
by the covalent crystal field may be detrimental to the 
calculation of the isotropic magnetic coupling. 

In fact, the titanates are notorious for the difficul- 
ties one encounters when trying to microscopically derive 
their properties. For example, Ref. Q finds a predomi- 
nant A-type antiferromagnetic coupling for LaTiOs while 
Ref. [3 predicts a ferromagnetic one, both contradicting 
the experimentally detected G-type coupling of that ma- 
terial. Our work on that compound [9j has yielded the 
correct magnetic order, but the application of the same 
model to YTiC>3 turns out to be not so successful. Similar 
problems have been reported in other studies of YTiC>3. 
Reference while deriving ferromagnetic couplings, 
predicts (in contradiction to the experiment) a strong 
anisotropy between the intra and the inter-plane cou- 
plings, i.e., J12 = —2.0 meV and J13 = —0.6 meV. Ref- 



erence [j| finds antiferromagnetic values for both these 
couplings in a parameter range which is considered to be 
realistic. Both these papers use models which are differ- 
ent than ours, but they also employ perturbation theory 
to second order in the Ti-Ti hopping to derive the re- 
quired superexchange parameters. 

It should be emphasized, however, that the starting 
point of our model, i.e., the crystal field and the or- 
bital ordering it implies do give a faithful description for 
YTiC>3. The failure of our model in producing correctly 
the isotropic Heisenberg coupling between the ab planes 
is likely to be related to the use of low order perturbation 
theory and to subtle inaccuracies in the parameters used. 
The alternative possibility suggested in Ref. Q based on 
orbital fluctuations is, in our opinion, not adequate, since 
it defies the experimentally-detected orbital order of the 
ground state. 

In view of the above difficulties, and since it is known 
that perturbation theory may be insufficient for the 
leading isotropic couplings but may well be reliable 
for the anisotropic ones, we have combined together 
the experimentally-deduced isotropic couplings of YTiC>3 
with the computed anisotropic ones, to calculate the clas- 
sical magnetic ground state. The result turns out to be 
satisfactory, when compared with experiment. Similarly 
to LaTiOs, we obtain a G-type moment along the crys- 
tallographic a axis, an A-type moment along the b axis, 
and a ferromagnetic moment along the c axis, the latter 
being the predominant one. Remarkably, this detailed 
structure is caused by the anisotropies, and cannot be 
derived solely on the basis of symmetry arguments. 

An even further check of our procedure is provided 
by the calculation of the spin-wave excitations. We find 
four dispersions: Three of them have considerable zone- 
center gaps, 6 meV, 12 meV, and 17 meV, while the 
fourth one has a very small gap, of the order of 0.3 
meV, and is approximately isotropic over the magnetic 
Brillouin zone. We have demonstrated that all branches 
have experimental counterparts as can be deduced from 
neutron-scattering data. Comparing the calculated dis- 
persion with the experimental one, we have found that 
they are in a plausible agreement. 
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